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ABSTRACT 

We present an inverse mapping approach to determining the emission height of the optical 
photons from pulsars, which is directly constrained by empirical data. The model discussed 
is for the case of the Crab pulsar. Our method, using the optical Stokes parameters, deter- 
mines the most likely geometry for emission including magnetic field inclination angle (a), 
observers line of sight angle (x) and emission height. We discuss the computational imple- 
mentation of the approach, along with any physical assumptions made. We find that the most 
likely emission altitude is at 20% of the light cylinder radius above the stellar surface, in 
the open field region. We also present a general treatment of the expected polarisation from 
synchrotron source with a truncated power law spectrum of particles. 

Key words: Isolated Neutron Stars, polarisation, Crab pulsar, inverse mapping, synchrotron. 



1 INTRODUCTION 

Despite over forty years of theoretical and observational studies, we 
still do not fully understand the emission mechanism responsible 
for the observed pulsar radiation. 

In an attempt to fully understand the high energy nonthermal 
emission from pulsars, both direct and indirect methods have been 
used. Direct models attempt to describe fully the source physics 
that subsequently generates the radiation we see. Conversely, in- 
verse models use the observed radiation in an attempt to infer in- 
formation about the source physics and emission geometry itself. 
In this paper, we present the methodology for a generic inverse 
method, whose aim is to restrict pulsar geometric parameters (a, \, 
and emission location), using few assumptions and many observa- 
tional constraints. Inverse models have had a long history within 
pulsar astrophysics. The models effectively utilise the link between 
lightcurve morphology and local source geometry, to place restric- 
tions on various source parameters, i.e., to restrict the pulsar's mag- 
netic inclination (a), the observer's viewing angle (x), and/or to 
determine the emission location within the pulsar magnetosphere. 
Perhaps the be st known of these models is t he Rotating Vector 
Model (RVM) (Radh akrishnan & Cookd [[969)), 

which constrains 

pulsar inclination and viewing angle. It does so by assuming that 
the radio linear polarisation is described by a polarisation vector 
fixed to local B field lines, (located close to the magnetic pole), 
which sweeps by an observer's line of sight as the pulsar rotates. 

Inverse methods in general, have evolved over time. Early 
ideas were to associate the width of pulse profiles with the rela- 
tivistically beamed opening angle of isotropic radiation, fro m point 
sources corotating within the pulsars magnetosphere (e.g.. ISmithl 
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1970), Zheleznvakov dl971 ), Izheleznvakov & Shaposhnikov 
1 19721) ). iMorinil d 19831) proposed a model for pulsed optical — s> 
7-ray emission based on relativistic electron beams producing 
highly beamed emission in the far (measured radially) magne- 
tosphere. Assuming emission from only the (purely dipolar) last 
open field lines of the equator of an orthogonal rotator, the ar- 
rival phase of emitted photons (assumed tangential to the local B, 
subject to aberration) were mapped onto observer phase. Align- 
ing this radius-to-phase mapping with the relative peak arrival 
phases for the Vela pulsar, Morini estimated that the radio, op- 
tical and 7-ray emission for the Vela pulsar, originated at dis- 
tances of ~ 0,0.5 and 0.7 Rlc from the neutron star surface. 
The advent of readily available computer power allowed the con- 
struction of radius-to-phase maps, as the assumption of tangen- 
tial beaming means that the mapping is independent of the spe- 
cific emission process. In t he lat e 1980s, Smith and co-workers 
dSmitfJ dl986h . lsmithetalj dl988l) ) extended this idea somewhat, 
producing radius-to-phase maps for the whole equatorial plane 
of an orthogonal retarded dipole. RY95 dRomani & Y adigaroglu 
( 1995)) succeeded in reproducing peak separations of the gamma 
ray pulsars, together with the radio-peak to gamma-ray peak rel- 
ative phase offset, while simultaneously restricting pulsa r inclina- 
tion and viewing ang le. iDvks. Rudak & Harding] d2004l) describe 
how light flight time delays can be respons ible for asymmetries in 
pulse profiles. DR03 dDvks & Ru dak 2003) developed a 'Two Pole 
Caustic' model (TPC), which can produce double peaked (and to 
a more limited extent, single peak) light curves. While certain re- 
strictive assumptions are sometimes necessary and justifiable, we 
have designed a general method tailored specifically at restricting 
unconstrained pulsar parameters. We effectively carry out a large 
scale computational search through pulsar parameter space to con- 
strain pulsar parameters. The aim was to create a methodology of 
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testing and refining hypotheses of how pulsars work. In this paper, 
we apply the model to pulsed optical emission from isolated neu- 
tron stars. 

Following a successful two years of operation, the Fermi 
Gam ma-ray observatory ha s identified over fifty gamm a-ray pul- 



Physical Component 



titty ) 

sars jAbdo etal.1 feOlObh & ISaz Parkinson et alj d2O10h . of which 
21 are radio-quie t gamma-ray pulsa r s and surprisingly, 11 are mil- 
lisecond pulsars ( Abd o et all d2009l) . iRansom et al. (2010)). From 



thes e observations, the polar-cap model has been effectively ruled 
out dAbdo et all d2009h . I Venter? Harding & Guillemotl d2009h ). al- 
though the polar cap will play a role in the initial particle accel- 
eration. This leaves the slot-gap, outer-gap (extended outer gap 
for millisecond pulsars ) and v a riations of the striped pu l sar win d 
dArons & ScharlemanrJ dl979h. ICheng. Ho & Rudermanl dl986id) . 
IPetril d2009h & lKirk, Skiraasen & Gallantl ( |2002|) ). as the main con- 
tenders for a theory of high-energy pulsar emission. 

It has long been recognised that the polarisation of pulsar ra- 
diation gives an indicat i on of the geometry of the emission zone. 
iThomas &G angadharal d2010h showed recently that the core and 
conal components of the radio emission from PSRs B 1839+09, 
B1916+14 and B21 11+46 were at low altitud es < 5% of the light 
cylinder radius. [Weltevrede & Wright! d2009h mapped the magne- 
topshere of PSR B 1055-52, indicating that the emission height was 
arou nd ~700 km above the neutron star's su rface. In a more general 
way, [Thomas, Gupta & Gangadharal f2oToh have looked the geom- 
etry for radio emission assuming curvature radiation. Our objective 
in this paper, is to investigate the use of optical polarisation data 
to map the regions within the magnetosphere at which the pulsar 
emission can originate. The advantage of optical radiation is that 
its emission mechanism, incoherent synchrotron, has a simple rela- 
tionship between its polarisation profile and the underlying geom- 
etry. It is also the only region where we have decent high-energy 
polarisation data, and such data is highly sensitive to the local emis- 
sion physics conditions. 



2 DESCRIPTION OF THE MODEL 

Historically, inverse models have focused on restricting various ge- 
ometrical parameters, such as pulsar inclination angle, viewing an- 
gle and/or the specific emission site responsible for the observed 
emission. These variables are assumed unknown initially and de- 
fine a parameter search space, where the aim is to restrict them by 
comparison with observation. However, all such models use rather 
restrictive initial conditions to constrain these parameters, for ex- 
ample, tangential beaming from specific subsections of the magne- 
tosphere using predefined emissivity functions. By removing these 
restrictive initial conditions and using modern compute power, we 
believe that the domain of applicability of inverse models can be 
increased. For example, we can apply the inverse approach to con- 
strain and test emission properties, as well as geometric properties 
of pulsar high energy emission. Effectively, one considers the vari- 
ous emission parameters as being undefined and adds these parame- 
ters to the overall search space to be constrained. Extending inverse 
modelling in this way may provide a new way of constraining pul- 
sar variables, and of approaching the pulsar problem in general. We 
describe in this section the overall approach we have designed and 
the specific inverse method we have developed to investigate pulsed 
nonthermal optical emission from pulsars. 

Our approach is composed of a number of conceptually dis- 
tinct steps, outlined in the flowchart in Figure [T] The main steps 
are subdivided into physical (P), computational (C), statistical (S), 



Assume synchrotron radiation from a 
power law distribution of particles 
within a dipolar pulsar magnetosphere 
[considering appropriate relativistic 
effects and light flight time] 



Computational Component: 



For discrete a with 0<a<180: 



/ "Mapping" from magnetosphere 
location to observer (a, 3C) 
polarised lightcurves ] 



Model magnetosphere as a grid of 
discrete points. At each point (r) simulate 
the above 'Physical Component' to create 
phase resolved polarised lighcurves: 

I (a, x, <J>) 1 

Q(a, X. *) 

U( a > X. *) 

y(a,x, *) J 

Also, record each grid point's polarisation: 
£(R. ct.z.*) 



= |(a,x,4>) 



5_ 



Statistical Component: 



[ Selection of which simulated 
lightcurves best fit the 
observations - thereby restricting 
(a, 30 ; 



Mapping / Inverse Mapping Component: 



I "Inverse Mapping" from 
phase resolved simulated 
polarised lightcurves to 
magnetosphere locations ] 



There is now a set of lightcurves, §(a,x,4>) 
for o <a <180, Carry out 
chi-squared tests on these lightcurves 
to determine which is closest to the 
observational data, 
[for a given pulsar, this should yield a 
'best fitting' (a A ,x B )combination] 



For the best fitting polarised lightcurves, 
I(«a.Zb.*) use |(K. ot A , x B .*) to 
to locate the emission region wherefrom 
the best fitting lightcurves originate. 



Figure 1. Representation of 'Inverse Mapping'/'Search Algorithm' ap- 
proach in flow diagram form. The approach is conceptually divided into 
physical, computational, statistical and mapping/inverse mapping compo- 
nents. 



and mapping (M) elements. Briefly, these steps constitute: (1) the 
application of a plausible physical model responsible for the emis- 
sion, (2) the computational implementation of this model, resulting 
in the creation of phase resolved lightcurves (£($, a, x)lE (3) the 
statistical comparison of the simulated lightcurves to observations, 
which results in a restriction of (a, x) parameter space enabling, 
(4) an 'inverse mapping' into the magnetosphere, locating the re- 
gions of emission responsible the best fitting lightcurves in the first 
place. 

Phase resolved Stokes parameters are produced for all 
possible values of 0° «S a sC 180° and 0° sC % < 180°, in finite 
increments Aa and Ax, based on the assumption that a simple 
radiative emission process is occurring globally within the mag- 
netosphere. The computational model divides the magnetosphere 
into a grid of points and the phase resolved angular distribution 

1 Emission from a given pulsar is simulated and phase resolved (<E>) 
lightcurves produced for a range of viewing angles (\). This process is re- 
peated using different inclinations (a) for the magnetic axis (which is not 
in general constrained from observations), f can represent numerous obser- 
vational properties, in this case the set of Stokes parameters I,Q,U,V. 
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of photons from each point (R) is recorded as £(R, $,a,x)- 
The summation of £ from all grid points creates a mapping from 
magnetospheric location to phase resolved observer profiles (as a 
function of (a, x))- The statistical component selects the £(a,x) 
which best fits the observations and a, x) allows an 'in- 

verse' mapping of the photon paths back into the magnetosphere. 

These steps constitute what we refer to as the 'Inverse 
mapping' or search algorithm approach to isolating the regions 
within the magnetosphere which may be responsible for pulsed 
optical emission. The main steps in the process are as follows: 

P-I Particle motion (outside of accelerated regions) is dom- 
inated by the magnetic field, making the field an important 
element in any model => we propose a magnetic field struc- 
ture in the form of a retarded d ipole, as the standard Deutsch 
formalism. lMichel & Lil dl999h . 

P-II It is plausible that synchrotron radiation is responsi- 
ble for pulsed optical emission, with the observed power law 
spectral form indicating that the underlying particle spectrum 
is also of this form. Using the premise of first order ap- 
proximations, together with no a-priori assumptions regard- 
ing favourable emission locations, we assume that each point 
in the magnetosphere can emit synchrotron radiation from an 
underlying particle power law, where the total number density 
of particles at any one point is equivalent to the Goldreich- 
Julian density (naj oc • B z (r)). 

P-III Each point R is allowed to radiate synchrotron radia- 
tion, which is characterised by the properties of the particle 
population at that point - essentially the energy distribution 
of the particles (a power law from P-II above), and the pitch 
angle distribution (hereafter PAD) of the particles. The pitch 
angle distribution is perhaps the least well known parameter 
of pulsar emission models. Since the PAD is not well con- 
strained, we assume various forms for the PAD and analyse 
the resultant lightcurves to choose which PADs may best rep- 
resent the emission. For all PADs, we assume that the emit- 
ting particles are symmetrically distributed about the mag- 
netic axis, with a cutoff occurring at a specific pitch angle, 
beyond which no emission occurs, PAD co . 

C-I Computationally, we represent the 3D volume of the 
magnetosphere as a grid. Each point on the grid is specified 
by local parameters such as (R, B, naj, v co ), where R, B, 
tigj, and v co are the spatial coordinate, the magnetic field, 
the local particle density and the corotational velocity respec- 
tively, at that point. The grid itself, is based on a spherical 
geometry, consisting of a series of fixed concentric spheres. 
The origin of the spheres is located at the centre of the neu- 
tron star, with points located such that R — (R, 0, (j>) in usual 
spherical coordinates. 

C-II At each grid point (R), a power law distribution of par- 
ticles, having a specific PAD, emits synchrotron radiation into 
a certain section of sky of solid angle d£l. Relativistic aber- 
ration, beaming and doppler shifts will alter the direction of 
emission, the irradiated solid angle and frequency of radia- 
tion, as seen by a sky based observer compared to an observer 
in the particle rest frame. These effects are considered and 
the emission from each point R, as seen by a sky based ob- 
server, is recorded in increments of Ax = 1°, for the range 
of 0° «S x < 180°. 

C-III Since one aim of this approach is to restrict both the 



inclination angle a and the viewing angle x, step C-II is re- 
peated for all inclination angles with 0° ^ a ^ 180°, in 
discrete increments of Act. 

S-I Steps C-II and C-III create a phase space (a,x), 
where each point in this phase space contains phase resolved 
lightcurves of the form: 



(i) 



with 0° < x < 180°, 0° s: a ^ 180°, and where £ rep- 
resents each of the four Stokes parameters (I ,Q,U,V). A 
X-squared goodness of fit test is performed between each of 
the simulated Stokes parameters and the observed Stokes pa- 
rameters, at each point in (a, \) phase space, to produce the 
best fitting (a, x) combination. 

M-I Knowing the best fitting (a, \) combination, 
£(R, 0, <f>, <£>, a, x) is used to inverse map the constituent pho- 
tons back into the magnetosphere, thereby locating the origin 
of the emission, which creates the best fitting lightcurves (£). 
The visualisation of these regions will hopefully yield an 
indication of the originating locations. 

In the following sections, each of these steps in illustrated in 
more detail, applying the search algorithm to a simulated pulsar 
with Crab like parameters: (P = 33 ms, P = 4.209 x 10" 13 ss _1 
and B SUI{ = 3.8 x 10 12 G). 



3 DESCRIPTION OF THE PHYSICAL MODEL 

It is assumed in this modelling that at any point in the magneto- 
sphere, radiation arises through synchrotron radiation from a power 
law spectrum of particles: the particle index, p, is fixed via the 
observed photon spectral index s, through the relation, s=(p+l)/2, 
where: 



F v oc v 



N(E) oc E- p . 



(2) 



When simulating emission over a range of frequencies (e. g., 
for U, B, V bands), contributions are summed only from those par- 
ticles whose energies, are such that they contribute significantly to 
emission at the observed frequency v. 

3.1 Synchrotron Radiation 

Synchrotron radiation from a single particle, spiralling around a lo- 
cal magnetic field (£?), at a pitch angle a (where n ■ B — cos a 
with h the instantaneous velocity of the particle), produces a con- 
tinuum spectrum of frequencies with a peak emissivity close to a 
critical frequency, f c where: 



fc(E,B, 



3 f ■ a( E 
2 fBSme {E- 



(3) 



with fs = eB/2irm, being the particle cyclotron frequency. 
To first order, the spectrum exhibits a quadratic rise and exponen- 
tial drop as the frequency passes through the critical value. At lower 
frequencies, the spectrum consists of discrete harmonics of the fun- 
damental (equation frequency below which no emission is pos- 
sible. 



ff 



Sb Eq 
sin 2 E ' 



(4) 
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Figure 2. Important source and observer geometric parameters relating to 
synchrotron emission. 



The emitted radiation is in general, elliptically polarised with 
the principle axes of the polarisation ellipse aligned parallel and 
perpendicular to the projection of B on the plane transverse to the 
emission direction, I (see Figure[2j. The angular dependence of the 
polarised emissivity depends strongly on the angle between the line 
of sight (I) and n, where we write, h ■ I = cos ip with sgn(ip) = 
a — 9. The case ip — 0, results in linear polarisation perpendicular 
to B, ip > denotes elliptical polarisation having the major axis 
perpendicular to the projection of B, and ip < indicates right 
handed polarisation and ip > left handed. 

For a population of particles, one must integrate emission 
over the velocity distribution, which amounts to an integration 
over the pitch angle distribution and over the range of particle 
energies present. Results for a power law spectrum of particles 
(N(E) oc E ') are well known; the emitted radiation spectrum 
is also a power law, with the photon spectral index related to the 
particle index, where a = (7 + l)/2 and 7 < 1/3, in order to 
maintain a finite integral. This monotonic spectrum is of course a 
special case, which results from integrating over a population hav- 
ing infinitely wide energy bounds. So that practically, the formal- 
ism can be applied only at frequencies which are unaffected by any 
issues related to the nature of finite particle energy bounds. 



3.2 A Truncated Power Law Spectrum of Particles 

In any real situation, one is dealing with particles having finite 
upper and lower energy bounds, which alter the spectrum signif- 
icantly from the classical monotonic scenario. For the purposes of 
this work, the angle dependant and integrated polarisation prop- 
erties of synchrotron emission from a truncated power law spec- 
trum of particles, with N(E) oc E~ J is considered. In this spec- 
trum, any given particle's energy is confined between lower (Ei) 
and upper (E2) limits, such that Ei ^ E ^ E2 having a (po- 
tentially) isotropic, axially symmetric (about B) pitch angle dis- 
tribution. Such a distributio n is rather general and has bee n dealt 
with in detail by GLW74 dGleeson. Legg & Westfoldll 19741) . They 
derive emissivity properties using the polarisation tensor for light, 
which represents the cross correlated quadratic components of the 
electromagnetic wave field as a rank 2 tensor, p a p (where a and 
/3 represent the component directions in the transverse plane). The 
tensor p a p, contains the complete polarisation content of the ra- 
diation and can therefore be related to the more commonly used 



Stokes parameters, as illustrated in equationQJo] For reference pur- 
poses, the underlying polarised emissivity properties are given here 
as equations [5] reproduced from GLW74. 



1 = ^(~r /2 ^)(^) C7+1)/2 



-(7-l)/2 



[7 = 



V = ^£(|)^ 2 0(e)cotc?(^ H sinc?)^ 2+1 i/- ( ^ 2) 
V 3 2 



ft {7 /2+i) + (1 + g(0))(C l/2 - -J l/2 ) 



(5) 



where: 



Mx) = / ?- 2 F(Qdt, Cn{x) = / C- 2 F 9 (QdZ 



(6) 



K n {x) = I C~'F*(t)d£, HO = J K 5/3 (y)dy (7) 



F P (0 = £K 2/3 (0, F„(0= ^i/3(0 



so that: 



{J TL 



C~ 3 F(0<% etc. 



(8) 



(9) 



where the K a , are modified Bessel functions and the variable x = 
f/fc where f c is the critical frequency given by [3] 

The effect of the truncated power law energy spectrum is 
seen in the fact that the functions modifying the observed Stokes 
parameters (i.e., J n (x), C„{x) and lZ„(x)), must be evaluated be- 
tween upper (a;i) and lower (£2) limits of the parameter 'x', where 
x = x(B, 0, E). Each polarised emissivity parameter is dependant 
on a functional combination of (f)(6), fs sin 9, f, 7, and this func- 
tionality is modulated by finite integrals over the particle population 
(having integrand 'x'). 

As analysed in detail by GLW, having a finite particle 
energy range results in an emission spectrum with lower (/;) 
and upper (f u ) frequency bounds. This spectrum possesses 
roughly two spectral breaks at frequencies f a and where 
fio < fa < /ft*' < fhi, with differing spectral indices for each 
polarisation parameter, each asymptotically dependent on the 
power law index 7 (see for example, Figure[3}. 

The frequency bounds, the internal transition frequencies and 
the form of the spectral variation, are all the result of energy cut- 
off effects. All effects have a simple physical basis and can be de- 
scribed in terms of the very useful variable of integration 'x', which 
contains the frequency of interest, f, and encodes all relevant phys- 
ical factors relating to the emission process, (E, B, 6). In fact, the 
fiducial frequencies are derived from the corresponding fiducial 'x- 
factors', designated xi < x a < x^ < x u . 

Physically, both the low and high energy frequency cutoffs 
are due to the absence of particles with energies, E, greater than 
the upper limit E2 (fio oc E^ 1 , fc oc En 1 )- Above x^\ we have 
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<B <B 



>E 2 




>E : 



\ = (E^)f, 



logf 

C = f B sm0(§K(y) 



Figure 3. Spectral variation of emission (Stokes I) from a truncated power 
law energy spectrum. The effect that decreasing the local magnetic field 
strength (B) and increasing the upper power law cutoff energy (E2) have 
on the spectral break points is indicated by the directional dotted arrows at 
the appropriate points. Our spectral range covers UBV. 



x'i <C x <C X2, so that the upper and lower limits of integration 
can be approximated to 00 and respectively. This corresponds to 
the classical scenario in which f c \ <C / <C f C 2\ the case in which 
the integrals are dependant solely on the spectral index 7, and can 
be expressed in terms of a combination of Euler Gamma functions. 
The photon indices for the stokes parameters are now the classi- 
cal values, the frequency dependencies already given in relations[5] 
i.e., a = (7 — l)/2 for Stokes I and Q, and a — 7/2 for Stokes 
V. Below x b , the finite limits of x\ and x-i cannot be ignored, and 
one gets a 'falling away' from the classical photon index as the fre- 
quency decreases. This results in an explicit dependence on x\ and 
X2 (and subsequently an altered frequency dependence). Deriving 
Power Law Expansion (PLE) approximations for the relevant in- 
tegral evaluations, GLW74 showed that the asymptotic spectrum 
should now vary as (f/f B sin0) 1/3 -E 7 ~ 1/3 for Stokes I and Q, 
and as E^ 1 for Stokes V (where E = E\ or E2). 

At fff \ again, the finite limits of xi and X2 cannot be ig- 
nored and one gets a 'falling away' from the classical photon index 
as frequency decreases. Power law approximations to the integrals 
in relations [5] show an asymptotic dependence onii, providing a 
spectral variation of f 1 ' 3 . 

As lower energy particles will have a higher fundamental 
frequency than higher energy particles, a point may come where 
some of the lower energy particles in the population cannot con- 
tribute to the frequency of interest, and the frequency at which 
this occurs, is designated f a . GLW74 deal with this case by de- 
termining the lowest energy particles from the population which 
can contribute to emission at frequency f (designated by E x where 
E 1 > _Ei), and replacing the upper limit of integration with x\, 
where x 1 — f/f c (Ei). It is this replacement, which results in the 
augmented spectral index of a — —7 at frequencies below f a . 

The different energy cut-off effects delineate different regions 
in which different physical effects alter the emission characteris- 
tics. The effects are encoded in the integration x-factors, both in the 
form of xi and the extent of the range (a;;, x%), over which emis- 
sion is allowed. Figure[3]shows the expected form of the emissivity 



for a truncated power law particle population emitting synchrotron 
radiation[f] 

In our approach to modelling emission from a truncated power 
law particle population, we follow the GLW74 approach and de- 
scribe the spectral variability in terms of correction factors, C\ 
i — 1,2,4 applied to the standard power law dependences (see 
AppendixfA}. We note that in limited circumstances, it is possible 
for x a > x^ , so that the spectral emissivity has a lower and up- 
per bound with a single internal transition frequency. This situation 
would arise when the upper and lower bounds are sufficiently close 
together. For completeness, we include the relevant correction fac- 
tors in AppendixfA](as this situation was not discussed in GLW74). 



3.3 Effect of source motion 

The description of the polarisation parameters above is valid when 
the magnetic field is stationary relative to an observer. In the case of 
emission from synchrotron radiating particles constrained to move 
along magnetic field lines in a pulsar magnetosphere, the underly- 
ing field structure is (relativistically) rotating. To describe the effect 
of the source motion on the observed polarised emission, we con- 
sider how a general Lorentz boost transforms the radiation fields. 
Also, to extract angular dependences and for frame transforma- 
tions, it is convenient to use the polarisation tensor, p a p, repre- 
sentation of light. 

In our notation, a noninertial observer in the corotating frame 
S (basis X yz), sees a stationary magnetic field B and views this 
field along a direction n, such that B / ■ n oc cos The observer 
sees light from a particle, if n is sufficiently close to i"(t), where r 
describes the trajectory of the charged particle (see FigurefJJ. The 
electromagnetic field of the light oscillates in a plane transverse 
to n (the 'observer plane' K), and may be expressed in component 
form along two mutually perpendicular axes (ii , 12) in the plane K. 
Frame S' observes the magnetic field B f to move at an arbitrary 
constant velocity (3 = v/c. An observer in S' will therefore see 
light emitted at an aberrated direction n' about the boosted field 
direction B' . The components of the electromagnetic wave will 
also be boosted (E' , B'). A general Lorentz boost, where frame S' 
moves at an arbitrary constant velocity (3 to frame S, is given as 
(equation 1 10b where a four-vector x transforms as a;" 1 = h^ L v x v . 



7 



-to 



-7A, 
-70. 



1 + 

"■- 
a' 2 



si Ac 2 



-to 



P 
1 + 



-to 

i + ^to 



(10) 

The electric and magnetic fields, transform as the components 
of the antisymmetric tensor F^u, (equation II lb. so that F'it = 
A. l3 k A a i F afi . The transformed wave four-vector (fe'), electric (E') 
and magnetic (B') fields are given in equations 1 1 2 1 — [T4l 



2 We note that the photon indices quoted above are asymptotic values eval- 
uated in the limit of 'small', using power law expansion representations of 
the functions J n {x), C n {x) and 1Z n {x) derived by GLW74. The transi- 
tion frequencies and photon indices are therefore asymptotic and indicative 
of the actual spectrum, which should transit smoothly through the transition 
frequencies and deviate smoothly from the asymptotic spectral indices. 
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It is now necessary to describe the observed intensity and po- 
larisation in terms of quantities seen in the observers frame (S'). 
For this, the polarisation tensor representation (p a p) of light is most 
useful, given by: 



Pa/3 = -J—& a C.p, 



(15) 



where £ a and Sp are two orthogonal components of the electro- 
magnetic wave field. The polarisation tensor can also be expressed 
in terms of the more widely used Stokes parameters, as in equa- 
tion[T6l 
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(16) 



Since the fields of the wave remain transverse in any frame, 
it is apparent that the tensor p a p will remain two dimensional in 
any new frame. Also, because the transformation is real, the real 
and imaginary parts of the tensor transform independently. These 
properties allow us to boost and rotate the electric field components 
from plane K in frame S, to plane K' in frame S'. Before boosting 
Pa/3, any reference to B in the expression for E' is removed, using 
the substitution B = [nE] . Also, a right handed orthonormal basis 
set (denoted Oj 1 j 2 j 3 ) is defined in the observers (S') plane (K'). 
The projection of the pulsars rotation axis (SI) on the plane K' is 
chosen as a uniform reference direction j\, so that the basis set 
designated as Oj 1 j 2 j 3 , can be defined through equation[771 



Ji 



[n'fl] 
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32 = [n Ji 
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(17) 



jt and j 2 are orthogonal and in the observers plane (K'), per- 
pendicular to the boosted emission direction j'3. We now express 
E' = Bi'/i + E 2 'j 2 , where Ex' = {E'j\) and E 2 ' = (E'j%). 
Since E = iE\i\ + E 2 i 2 , it is possible to write E' in the Oj^ 
basis explicitly in terms of Ei and E 2 as follows: 



(EiJitKji) = {{iEixn + E 2 X2i}j 1 , 
{iE 1 x 12 + E 2 X22}j 2 ), 



(18) 



where the x factors (xn etc.), are given by equation [T9j and 
are dot products between the Oj 1 j 2 basis vectors and the E\ and 
E 2 components of E' . 
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Equation [T8] expresses how the components of E and B in 
plane K, written in terms of the basis (ii , i2), transform into com- 
ponents of E' and B' in plane K', expressed along a basis j 1 ,j 2 - 
Using equation [T6] the components of the emission-polarisation 
tensor p aB in frame S" can be determined as follows: 



pu — £iS 1 xnx 22 + £ 2 £ 



2 X21X 2 1, 
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(20) 
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The polarisation tensor can be related to the Stokes parameters 
in the standard way (see equation[T6t. and using these relations, the 
Stokes parameters observed in frame S', can be expressed in terms 
of the field components of (E, B), as calculated in frame S. In the 
same way, (/', Q' , U' , V) in S' can be related to (7, Q, U, V) in 
frame S. Carrying out the algebra the results are as follows: 



V = 



2 {(^1 + * 2 i2){i + Q) + {'4i + xL)V - Q)} 

\ {{xlx - x\ 2 ){I + Q) - {x\ 2 - x 2 21 )(I - Q)} 



(21) 



{xi2Xu)(I + Q) + (x 22 x 2 i){I - Q) 

V = (x 22 Xn - X 12 X2l) V 

A significant source motion will alter the observed polarisa- 
tion, the intensity and magnitude of the Stokes parameters will 
change, and this will be accompanied by a Doppler shift in the ob- 
served frequency. The appearance of a non-zero Stokes U in the 
observer frame is not in itself significant, as the relative magnitude 
of Stokes Q and U, depend on the orientation of your reference axes 
in the observer plane. Figure [4] shows the source and observer po- 
larisation spectrum for a source magnetic field of 1.9 x 10 G, with 
a source velocity of j3 = 0.87. 



4 DESCRIPTION OF THE COMPUTATIONAL MODEL 

The above physical emission model is simulated within the pul- 
sar's magnetospheric environment. The three dimensional magne- 
tosphere is represented computationally in spherical polar coordi- 
nates, and it is assumed that emission arises from within the open 
volume of the magnetosphere, which is determined numerically. As 
such, emission is only simulated from those grid points lying in the 
open volume, where the grid itself is specified using coordinates 
(r, 8, <j)). The number of grid points on any given spherical surface 
is controlled by the parameters Ne and N$, (the number of 8 and 
(f> divisions). The parameter N r specifies the number of divisions 
in the radial direction, and the grid itself is devised in such a way 
that each grid point represents an equal volume of space, so that in 



Pulsars 7 




REL I: Stakes I 

REL Q: Stokes Q - 
REL U: Stokes U 
REL V: Stokes V 



■ 10 1e+11 1e+12 1e+13 !e+14 1e+!5 1e+16 1e+17 1e 



1e+09 1e+!D 1e+11 1e+12 1e+13 1e+14 1e+15 1e+16 1e+17 1e+18 



Figure 4. Effect of a relativistic source motion on the observed Stokes parameters, determined here for B = 1.9 X 10 G, /3 = 0.87. 



essence, the discretisation of the magnetosphere is homogeneous 
with respect to volume. 

The emission is calculated at all points in the open volume, 
which requires the determination of local parameters. These in- 
clude the magnetic field strength, B, the total number of particles 
at a specific point, naj, the assumed pitch angle distribution, lo- 
cal co-rotational velocity and the range of particle gamma factors 
which contribute to observable emission at that particular point. All 
necessary dependencies are calculated at each grid point to create 
the local emission profile, which effectively determines how the po- 
larised spectral emissivity varies across the extent of the local pitch 
angle distribution. Each profile is subject to the required frame 
transformations (described above), which determine the emissiv- 
ity profile as seen by a stationary observer external to the rotating 
magnetosphere. 

The pitch angle distribution at any open volume point, 
generates emission which extends over a finite solid angle. 
This solid angle of emission has a finite extent latitudinally 
and longitudinally with respect to an observer, which directly 
correspond to a range in observer viewing-angle, x an d temporal 
phase, respectively. Discretising this observer spac^f] allows 
emission from any point to be recorded as a function of viewing 
angle and phase. For any given model pulsar (i.e., for a specific 
magnetic inclination, a), the code records the phase resolved 
polarimetry in the form of an array, denoted as (,(R, 6, <f>, <&, a, \) 
(= £,(R, 4>, a, x))- This means that for each grid point location, 
R, the intensity and phase resolved Stokes parameters for the 
range of viewing angles over which emission is seen, are recorded^ 

The indices in the £ array are discrete and represents a record 
of the phase resolved Stokes parameters seen by observers from 
point R. Once all grid points have been sampled, the array £ can 
be summed from all appropriate locations to form ^ (ct, "^iS - 
i.e., the phase resolved Stokes parameters observed from the entire 
magnetosphere for a pulsar inclination of a. 

A single run of this code will create files which contain 



3 Representing 0° ^ \ ^ 180° and ^ $ ^ 1 with step-sizes of Ax 
and A$. 

4 The magnetosphere is of course rotating, so at any point in time, the 
correct phase of emission is determined taking light flight time, aberration 
and rotation into account. 

5 Represented in Figure[T]as the array £ (a. a i XB > 



£(a,x,&) f° r a single, specific pulsar inclination a. The search 
algorithm approach attempts to restrict a and x by simulating emis- 
sion for the range of all possible a valuefl (0° ^ a ^ 180°), and 
subsequently selecting those simulated lightcurves which best fit 
the observational data (which should correspond to a unique com- 
bination of a and x)- To generate lightcurves for different possible 
inclinations of a given pulsar, a is varied from 0° ^ ct ^ 180° 
in discrete steps Aa. This generates a large pool of phase resolved 
model lightcurves, each dependant on a unique combination of a 
and x - f° r example, with Aa = 10° and Ax = 1°, one gener- 
ates 180 x 18 = 3240 distinct lightcurves (x4 when polarimetry 
is considered). 



5 RESULTS 

The search algorithm approach was partly motivated on the 
premise that any simulated Stokes parameters (£), would be a 
variable function of the phase space (a, x, PAD, PAD co ). The 
simulated Stokes parameters vary as a smooth function of all 
phase space variables £(a, x, PAD, PAD C J3 ), thereby allowing 
for the selection of a best fitting parameter set. This being the 
case, one can restrict the pulsar geometry ((a,x) combination), 
which best fits observations and subsequently carry out the inverse 
mapping step to locate the associated source region within the 
pulsar magnetosphere. 

Simulated Stokes parameters are by definition dependent 
on all phase space variables (a, x, PAD,PAD co ). As described 
in section [4] this phase space is sampled discretely, thereby 
generating a large amount of model lightcurves. We find that the 
simulated Stokes parameters (£), vary smoothly as a function of 
this parameter space, showing systematic behaviour as we vary 
the observed energy band, viewing angle x, pulsar inclination a, 
and particle pitch angle distributions. We discuss initially, the form 
of the simulated Stokes parameters, how these vary as a function 
of parameter space and how this variation can be used to select a 
best fitting parameter combination (the 'search algorithm'). We 
subsequently describe our restrictions on pulsar geometry ((a, x) 



6 By explicit assumption, there is no a priori restriction on a (which is 
generally the case, as observations constrain neither a nor x). 

7 PAD refers to the pitch angle distribution which can either be isotropic 
or have a Gaussian structure; PAD co defines the maximum pitch angle. 
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Discretisation of observer's 
viewing angle and phase 



Figure 5. Centre: The concentric spherical Cartesian grid imposed on the pulsar magnetosphere, showing points only in the open volume of an aligned rotator 
(tl ■ fi = 1). Right: At each individual location (R), we calculate the radiation emitted from the local particle distribution, taking the local physical conditions 
into account - e.g., B, v co , the power law particle distribution, the specific pitch angle distribution etc. Left: Representation of how emission from a particle 
pitch angle distribution extends over a range of viewing angles (x) and phase (<J>). Emission is recorded computationally into discrete bins of extent A% X A$ 
(as described in section[4). 



combination), and discuss the inverse mapping step which lo- 
cates the associated source region within the pulsar magnetosphere. 

We note here, that our simulations resolve pulsar phase into 
200 bins, and we subsequently quote the phase extent of a single 
pulsar cycle as going from — > 1 . As choice of a reference phase 
is arbitrary, we define phase to be the phase at which the mag- 
netic pole associated with the quoted a value would be seen by an 
external observer. From symmetry, phase 0.5 is the phase at which 
the magnetic pole at a = 180° — a would be seen. We refer to the 
magnetic pole at inclination a as pole 1 and the magnetic pole at 
180° - a as pole 2. 



5.1 Form and variation of simulated Stokes parameters 

The search algorithm method, in sampling such a large parameter 
space, inherently creates a large amount of data. It is instructive 
initially to discuss the overall form of the simulated Stokes param- 
eters and to describe how they vary as a function of the phase space 
parameters a,x, PAD, and PAD co . Some useful systematic trends 
are evident which will be discussed. 



which should reflect the overall general trends in the variation of 
5.1.2 Stokes 1: Single and Double Peaks 

We have found that simulated lightcurves only have a single peak 
for any inclination a < 50°, but a double peak profile is seen to 
emerge for a > 50° . The intensity of the secondary peak increases 
as ot — > 90°, where the main: secondary peak ratio varies smoothly 
with viewing angle. This allows for an immediate restriction in 
terms of predicted inclination and viewing angle for single and 
double peaked pulsar profiles. 

Using £=1(R, a, \, it is found that the main emission 
contribution to the secondary peak, comes from pole 2 if \ ^ 90°, 
and from pole 1 if \ ?S 90°, indicating that as a > 50°, the 
secondary peak is formed from the magnetic pole opposite the pole 
contributing to the main peak. Figure|7]shows some phase resolved 
plots of the I(pole l):I(pole 2) ratios for an inclination of a = 80°, 
at different viewing angles, which illustrates the relative contri- 
butions of emission over each pole to the final integrated lightcurve. 



5.1.1 PAD and PAD co 

It is found that similar trends exist for the variations of £ between 
the different PAD combinations (i.e., each of the PAD functions 
based on isotropic, circular, linear and Gaussian profiles all have 
similar £(a, x))- We find that varying the PAD co has a larger effect 
on the polarisation than varying the PADitself. The general trend 
is for larger PAD co values to generate broader peaks, but the peaks 
are located at the same phase locations so that the overall variation 
of £(a, x) is similar for all PADs. Given this degenerate behaviour, 
we will focus our analysis on the isotropic pitch angle distribution, 



5.2 General Observational Restrictions 

Certain features of Crab pulsar emission are clear and distinct and 
any successful modelling of this object should be able to recreate 
at least some of these main features: 

[a] Double peaks: The Crab pulsar exhibits a double peaked 
structure at all wavelengths. Only simulated dipole inclinations 
with a > 50° are capable of producing double peak emission. 
This is our first restriction of phase space. 
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Figure 6. Relative intensity (Stokes I) at all viewing angles (V. A.) for successive pulsar inclinations. Parameter sets are (clockwise from top left) (a) = 
70°, 80°, 90° with (PAD, PAD co )=(isotropic, 20°) in each case. The secondary peak increases its intensity and viewing angle extent, as a — > 90°. At 
a = 90°, the 'secondary peak' is of equal intensity to the main peak. In terms of restricting parameter space, this evolution shows that the observed 
maimsecondary peak ratio of ~ 0.3 (for the Crab pulsar), can only be reproduced for a restricted range of viewing angles (for A% ~ 5° about a central x ) 
for a > 60° 



[b] Relative peak intensity and bridge emission: The double 
peaked structure is accompanied by the requirements of explaining 
the relative intensity of the peaks. A secondary peak whose 
intensity and extent of visibility (range of \ values over which it is 
visible) is a (monotonic) function of inclination a, appears within 
the simulations for a > 50°. The ratio of main: secondary peak in- 



tensity of 1:0.3, is possible for a narrow range of x associated with 
any a > 60°. This range of \ extends no more than ~ 2.5° about 
a central \ — X°- We find that \o is a smooth function of pulsar 
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Figure 7. Intensity observed for (a) = 80° showing relative relative flux 
contribution of emission from pole 1 and pole 2 for different viewing angles. 
Here pole 1 refers to the magnetospheric regions above pole 1 and not just 
to the polar cap regions. Results illustrate that, depending on whether \ $ 
90°, different poles will contribute to main and interpulse emission. 



inclination (\o = Xo (<*)), Xo decreasing steadily as a increases(f| 
Our second restriction can be summarised as follows: only for 
a > 60° and for viewing angles Ax ~ 2.5° about \ — Xo(a), 
are simulated lightcurves able to reproduce the maimsecondary 
peak ratio of ~ 1:0.3. Bridge emission is also a function of a 
where the large bridge emission seen at lower inclinations (a), 
tends to disappear for the more orthogonal rotators. 

[c] Peak phase separation: The Crab pulsar's peak separation is 
~ 0.4 in phase. Results indicate that the simulated model is unable 
to produce a double peak structure separated by less than 0.4 in 
phase, where in fact most viewing angles have peak separations 
of 0.5. In hindsight, this is a result of the inherent symmetry of 
the physical model, as it simulates emission from the open volume 
of both poles of a symmetric dipolar magnetic field structure. 
The influence of relativistic beaming can lead to asymmetries in 
pulse structures, so one would sensibly infer that a global emission 
model viewed at a 7^ 90° , x 7^ 90° should subsequently exhibit 
non-symmetric behaviour. It seems therefore, that the explanation 
for the symmetry in these model results, most likely lies in the 
averaging of emission over such a large volume. The large width of 



Xo > 80° for a = 60°, \o ~ 47° for a = 80°, as can be seen in 
Figure [6] 



the peaks is also most likely due to the large volume of emission. 
It is hoped that the inclusion of more general physical constraints 
into the model, should lead to a removal of the symmetry. 



5.3 Inverse Mapping - A two pole emitter? 

Our results given that the model results are capable of more suc- 
cessfully reproducing lightcurve flux profiles than the associated 
polarisation profiles, inverse mapping results are given for a se- 
lected (a, x) combination, which compares most favourably to the 
observed double peak structure of the Crab pulsar and the observed 
QU relationship. A value of (a, x) = (70°, 45°) is chosen for 
(PAD, PAD co )=(isotropic, 20°), these values were bas ed on a best 
fit in t he x 2 sense between the simulated data and the lSmith et al.l 
( 1988) observations. In Figure [8] we show the intensity and de- 
gree of linear and circular polarisation as a function of phase. We 
note that the inte nsity distribution is wider than normal observed 
- see for example Slow ikowska et al.l (2009) this can be explained 
by our initial approach of filling t he ope n field region with radiat- 
ing particles. Istowikowska et all d2009l) also showed an increase 
in the degree of linear polarisation happening prior to the main 
pulse which we see although not at exactly the same phase rela- 
tionship. From Figure|9]we see the expected swing in polarisation 
a ngle around the two peaks and QU plo t s mor phologically similar 
to lSlowikowska et al.l j2009h and lSmithl dl98ot). Our prefer r ed val - 
ues of a and x ^ broadly consistent with Ng & Romanil d2004h . 
which provided a robust estimate of the pulsar inclination angle. 

The inverse mapping element of the search algorithm allows 
us to decompose light-curves, with the option of selecting different 
phase resolved regions (on our current resolution phase is divided 
into 200 bins), which we can subsequently retrace back into 
the magnetosphere , isolating the originating locations of the 
component photons. This will localise the magnetospheric regions 
which contribute to emission at a particular phase, one of the main 
motivations for the search algorithm. 



We try and present the results of the inverse mapping in an 
interpretable sense here in 2 and 3 dimensions, although the task 
is not an easy one. The selected (a, x,PAD,PAD co ) combination 
chosen to illustrate the inverse mapping is (70°, 52°, isotropic, 
20°), which has a main peak at phase ~ 0.1 and an interpulse at 
phase ~ 0.6. We choose specific phase regions (A$ = 5) centred 
on the arrival phase of the main peak and the inter-pulse, so as to 
obtain a view of the magnetospheric regions which contribute to 
'peak' intensity. We note however on the basis of Figure [7] that the 
fit is a smoothly varying function of a, x,PADand PAD co . In Figure 
I5.3| we show the effect of changing the PAD co showing intensity vs 
phase for PAD co values of 1° and 5°. 

Carrying out the inverse mapping step for (a, \) = 
(70°, 45°), FigurelTTIpresents a 3D colour coded intensity map of 
the originating location of emitted radiation seen by an observer at 
X = 52° relative to the model pulsar. Figure QT] (left) shows the 
origin of emission which composes the main peak and Figure [TT] 
(right) shows the origin of emission composing the interpulse 
(as an aid to perspective, the circles define the light cylinder 
distance and the brownish field lines are the boundary of the open 
volume/closed volume region within the pulsar magnetosphere). 
Immediate points of note, are that the main peak is primarily 
composed of emission from one pole (which we term pole 1), 
while the interpulse is composed of emission from the opposite 
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Figure 8. Phase plots for 70° inclination and viewing angles 45° and 20° . We show the total optical intensity, and degree of linear/circular polarisation. 



pole (here defined as pole 2). Also, the most intense emission is 
located in tightly constrained volumes in the inner magnetosphere. 

It is immediately apparent that these results show localised 
emission regions and clearly indicate a two pole emission model. 
Figure QT] illustrates rather dramatically, that the first order as- 
sumption of a global truncated power law synchrotron emission 
model (even for a relatively wide pitch angle distribution of 20°) 
can produce a rather limited spatial extent of maximum emission 
within the pulsar magnetosphere, while matching certain features 
observed in the Crab pulsar profile. The main location of emission 
is well away from the standard outer gap location and closer to the 
polar cap centre, but a closer analysis (described subsequently) 
yields that the emission is at distances of ~ 0.2Rlc from the 
polar cap, further from the polar cap surface than polar cap models 
generally estimate. 

In an independent analysis, we attempt to localise relative 
emission contribution as a function of magnetospheric location. 
We subdivide the magnetospheric volume into concentric spherical 
shells of width Ap (p being the spherical radius), and also into 



coaxial cylindrical shells of thickness Ar (r in cylindrical coordi- 
nates). We integrate emission from all points within these shells to 
obtain a radial profile (spherically and cylindrically respectively) 
of absolute contributions to the total emission. The spherical 
profile is presented in the left hand panels of Figures IA1I and IA2I 
with the cylindrical contributions shown on the right. Figure IAT1 
shows emission contributions to the main emission peak at phase 
~ 0.1, Figure. lA2l shows emission contributions to the secondary 
peak at phase ~ 0.6. The top panels in these figures show the total 
emission from over both magnetospheric poles, the middle panels 
show emission regions associated with pole 10 only and the bottom 
panels show emission from pole j"1 only. 

From the previous figures, certain trends are apparent. Pole 1 
emission dominates for the main pulse and pole 2 emission for the 
secondary, each pole having distinct emissivity trends, both radially 

9 Pole 1 corresponds to the magnetospheric region above the magnetic pole 
at an inclination a 

10 Pole 2 refers to emission associated with the pole at inclination 180° + 
a. 
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Figure 9. Illustrative polarisation angle vs phase and QU plot for 70° inclination and viewing angles 45° and 20° 
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Figure 10. Phase and intensity relationship for 70° inclination and viewing angle 45° for two PAD of 1 and 5 showing the ragged nature of the intensity for 
the lower pitch angles. 



(cylindrically) and spherically. The poles dominating emission tend 
to have their maximum contributions located at very small spher- 
ical and cylindrical (~ 0.2Rlc) radii, whereas the poles which 
contribute less to the total emission are located further from the 
star. The high peak in the emissivity, may be due to the influence 
of higher local charge densities lower down in the magnetosphere. 
Generally, the collective interpretation is that the main contribution 



to emission is very localised closer to the neutron star surfacJH 
than to the light cylinder. This tends to agree with models attempt- 
ing to explain optical emission through a localised effect such as 
the slot gap. Inde ed our results are consistent with the emission alti- 
tude predicted bv lDvks. Wright. & Demorestl ( feOld) . Furthermore 
we note that the slot-gap models primarily localise emission in the 



11 The poles contributing less to the emission of a given peak also seem to 
have their emission localised, but at a greater distance from the star. 
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Figure 11. Three dimensional representation of the regions contributing to main peak (left) and secondary peak (right) emission for a model (Crab) pulsar 
with parameters (a, x,PAD,PAD co )=(80° , 25°, isotropic, 20°). The intensity scale shows the emission contributions at the regions mapped. For reference, 
the circles define the light cylinder boundary and the light brown lines are the open volume/closed volume boundary. Emission is seen to be concentrated in 
the lower magnetosphere from a single pole. 



transverse (co-latitude) direction whereas we also restrict emission 
in altitude. 



5.4 Conclusion 

Optical polarisation studies can be used to determine the local ge- 
ometry of the emission region. In particular from this work we have 

• A simple synchrotron model for the emission gives reasonable 
agreement with observations. 

• A prediction that the emission is low in the magnetosphere at 
an altitude in the region 30-40 stellar radii, that is away from both 
the polar cap and likely outer gap regions. 

• The linear polarisation peaks on the rising edge of th e 
main pulse, consistent w ith the findings of ISmith et al] dl988l) : 
ISlowikowska et alj d2009h . 

• The radiation for the most part is circularly polarised with an 
interplay between linear and circular polarisation on the rising edge 
of the main pulse. On this basis observations of pulsed circular po- 
larisation from optical pulsars would provide a significant geomet- 
rical restrict on pulsar parameters although we accept more detail 
work is needed of more realistic emission zones. 

• Pulse widths which are significantly greater than those ob- 
served although this probably stems from our requirement to fill 
to open field region with an emitting plasma. 

• Due the inherent symmetry of this model we do not see any 
significant bridge emission at the preferred orientation and viewing 
angle. 

Our future work will entail restricting emissions to localities 
around the last open field lines thereby removing the inherent sym- 
metry existing within the model. In this way we hope to able to 
make firmer predictions in relation to the optical emission specifi- 
cally and more generall y comment on correlations between the op- 
tical and radio emission ISIowikowskaet al. (2009). As well as the 



emission location and pulsar orientation our approach can also re- 
strict the PADand PAD co . We also intend to make more detailed 
comparison between these predictions and measurements of the 
pulsed circular polarisation. To-date there have been no measure- 
ments of optical circular polarisation from any pulsar. A new in- 
strument, the Gal way Astronomical Stokes Polarimeter (GASP), 
lKvneetalJl20Tcl) has this capability and observations are planned 
to observe the Crab pulsar in late 201 1. 
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where, 
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where x[ is the frequency limit cutoff x-factor as discussed in 
section 3.2. 

The correction factors encode the full details of the effects 
due to the truncated nature of the particle spectrum. They include 
the explicit dependency on xi and x 2 and further, they encode the 
physical effect of the frequency limit cutoff energy by replacing the 
limit x\ with, x\ whenever E\ < E[. 

Once x a < x^\ the four fiducial x-factors defining the end- 
points and transition-points are x t < x a < x\/ < x u . Using 
asymptotic approximations to the PLEs (Power Law Expansions) 
of GLW74, the variation of the correction factors can be expressed 
as follows, where i = 1, 2 and j = 4: 
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APPENDIX A: SYNCHROTRON POLARIMETRY FOR A 
TRUNCATED POWER LAW SPECTRUM OF PARTICLES 

The effect of the truncated power law energy spectrum re- 
sults in the functions modifying the observed Stokes parameters 
(i.e., J n (x), £„(x) and TZ„(x)), needing to be evaluated between 
upper (xi) and lower (x 2 ) limits of the parameter 'x', where x — 
f/fc- This introduces an explicit dependency on the underlying 
particle energy limits, Ei and E 2 . It is the functional form of these 
x-factor limits and their variation, which dictates the resultant po- 
larised emissivity spectrum for synchrotron radiation. 

As described by GLW74, the polarised emissivity can be in- 
vestigated by defining and introducing unique correction factors, 

(xi , x 2 ) for each of the Stokes parameters, where: 
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with ii ^ ij asii ^ ij. The correction factors have three major 
regions (referred to as Rl, R2 and R3 above) with asymptotically 
distinct behaviour. The boundaries are defined by the transition x- 
factors, which also define the frequency transition points in the pho- 
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ton spectrum. subsequently divides the spectrum into two regions as follows: 
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The photon spectrum corresponding to the transition x-factors 
are as follows: 
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S 4 representing Stokes 



There is one other possible spectral form which occurs if 
x* > x b z (7) and in this case the photon spectrum is composed 
of a single transition point, designated f a . One now has x-factor 
divisions, where xi < x a < x u - 

For xi < Xi < x a , part of the particle population is excluded 
from contributing radiation at the frequency of interest, because 
their fundamental emission frequency is too high. We label this re- 
gion Rib, as it corresponds directly to region Rl of the previous 
case (where x\ > a!j (7)). The range where x a < x\ < xm, we 
label region R3b- In this case the limits (xi, x2) are wide enough 
apart, so that we can relax to the classical scenario where the lim- 
its are approximated to (0, 00) respectively which corresponds to 
region R3 of the previous case. 

If x\ > x b l \-y), the spectral form of the correction factors 
and resultant emissivity in regions Rib and R3b, is identical to the 
spectral variation in region Rl and R3 shown in relations lA4l - |A7| 
The major effect of having x* > k*' (7), is to have a single inter- 
nal transition point (x a ) differing from both x a and x b (7), which 
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Figure Al. Relative integrated emission strength within the pulsar magnetosphere for a phase extent of 0.025, centred on the main peak (at phase ~ 0.05), 
for simulation parameters (a, \, PAD, PAD co ) = (80°, 25°, isotropic, 20°). Left panels: An integration of emission carried out within spherical shells 
concentric with the neutron star, each shell of thickness Ap = Right panels: Integration of emission using co-axial cylindrical shells of thickness 

Ar = R L c /50. Bottom panels show the location of emission associated with pole 2 only, middle panels from pole 1 only and top panels from both poles 
(see text for definition of 'poles' and a further explanation). 
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Figure A2. (As figure lATl but for secondary peak). Relative integrated emission strength within the pulsar magnetosphere for a phase extent of 0.025 centred on 
the secondary peak (at phase ~ 0.55), for simulation parameters (a, \, PAD, PAD co ) = (80°, 25°, isotropic, 20°). Left panels: An integration of emission 
carried out within spherical shells concentric with the neutron star, each shell of thickness Ap = Rlc /50. Right panels: Integration of emission using 
co-axial cylindrical shells of thickness Ar = Ft lq /50. 



